Set up environment.

Fraxinus americana

species <- 'fraxinus_americana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Amer'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus americana '),
                                                        'occurrences')), 
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus americana, zoomed out

species <- 'fraxinus_americana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Amer'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus americana '),
                                                        'occurrences')), 
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus caroliniana

species <- 'fraxinus_caroliniana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Caro'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus caroliniana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus caroliniana, zoomed out

species <- 'fraxinus_caroliniana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Caro'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus caroliniana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus cuspidata

** Check values in E Texas. plot(raster) extract(raster, click(n=1))

species <- 'fraxinus_cuspidata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Cusp'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus cuspidata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus cuspidata, zoomed out

species <- 'fraxinus_cuspidata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Cusp'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus cuspidata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus greggii

species <- 'fraxinus_greggii'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Greg'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus greggii '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus greggii, zoomed

species <- 'fraxinus_greggii'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Greg'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus greggii '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus nigra

species <- 'fraxinus_nigra'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Nigr'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus nigra '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer= ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus nigra, zoomed out

species <- 'fraxinus_nigra'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Nigr'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus nigra '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer= ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus pennsylvanica

species <- 'fraxinus_pennsylvanica'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Penn'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus pennsylvanica '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer= ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus pennsylvanica, zoomed out

species <- 'fraxinus_pennsylvanica'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Penn'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus pennsylvanica '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer= ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus profunda

species <- 'fraxinus_profunda'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Prof'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus profunda '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer= ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus profunda, zoomed out

species <- 'fraxinus_profunda'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Prof'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus profunda '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer= ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus quadrangulata

species <- 'fraxinus_quadrangulata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Quad'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus quadrangulata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer= ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus quadrangulata, zoomed out

species <- 'fraxinus_quadrangulata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Quad'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus quadrangulata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer= ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Now, for Lorenz ecbilt gcm.

Fraxinus americana

species <- 'fraxinus_americana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Amer'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus americana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus americana, zoomed out

species <- 'fraxinus_americana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Amer'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus americana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus caroliniana

species <- 'fraxinus_caroliniana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Caro'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus caroliniana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus caroliniana, zoomed out

species <- 'fraxinus_caroliniana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Caro'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus caroliniana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus cuspidata

species <- 'fraxinus_cuspidata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Cusp'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus cuspidata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus cuspidata, zoomed out

species <- 'fraxinus_cuspidata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Cusp'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus cuspidata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus greggii

species <- 'fraxinus_greggii'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Greg'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus greggii '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus greggii, zoomed out

species <- 'fraxinus_greggii'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Greg'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus greggii '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus nigra

species <- 'fraxinus_nigra'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Nigr'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus nigra '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus nigra, zoomed out

species <- 'fraxinus_nigra'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Nigr'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus nigra '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus pennsylvanica

species <- 'fraxinus_pennsylvanica'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Penn'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus pennsylvanica '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus pennsylvanica, zoomed out

species <- 'fraxinus_pennsylvanica'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Penn'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus pennsylvanica '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus profunda

species <- 'fraxinus_profunda'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Prof'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus profunda '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus profunda, zoomed out

species <- 'fraxinus_profunda'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Prof'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus profunda '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus quadrangulata

species <- 'fraxinus_quadrangulata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Quad'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus quadrangulata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus quadrangulata, zoomed out

species <- 'fraxinus_quadrangulata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Quad'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus quadrangulata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Now, for Lorenz ccsm gcm.

Fraxinus americana

species <- 'fraxinus_americana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Amer'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus americana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus americana, zoomed out

species <- 'fraxinus_americana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Amer'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus americana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus caroliniana

species <- 'fraxinus_caroliniana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Caro'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus caroliniana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)

Fraxinus caroliniana, zoomed out

species <- 'fraxinus_caroliniana'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Caro'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus caroliniana '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus cuspidata

species <- 'fraxinus_cuspidata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Cusp'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus cuspidata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus cuspidata, zoomed out

species <- 'fraxinus_cuspidata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Cusp'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus cuspidata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus greggii

species <- 'fraxinus_greggii'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Greg'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus greggii '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)

Fraxinus greggii, zoomed out

species <- 'fraxinus_greggii'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Greg'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus greggii '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus nigra

species <- 'fraxinus_nigra'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Nigr'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus nigra '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus nigra, zoomed out

species <- 'fraxinus_nigra'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Nigr'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus nigra '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus pennsylvanica

species <- 'fraxinus_pennsylvanica'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Penn'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus pennsylvanica '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)

Fraxinus pennsylvanica, zoomed out

species <- 'fraxinus_pennsylvanica'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Penn'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus pennsylvanica '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus profunda

species <- 'fraxinus_profunda'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Prof'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus profunda '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus profunda, zoomed out

species <- 'fraxinus_profunda'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Prof'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus profunda '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)

Fraxinus quadrangulata

species <- 'fraxinus_quadrangulata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Quad'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(rangeMap, border = 'blue', main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus quadrangulata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(envMap, add = TRUE)
plot(rangeMap, border = 'blue', add = TRUE)
map("state", add = TRUE)
map("world", add = TRUE)
points(records$longitude, records$latitude, pch = 16, cex = 0.6, col = 'red')

Fraxinus quadrangulata, zoomed out

species <- 'fraxinus_quadrangulata'
load(paste0('./cleaning_records/', species, '_buffer.rData'))
speciesAb_ <- 'Frax_Quad'
load(paste0('./Models/Maxent/model_outputs/', speciesAb_, '_GCM', gcm, 
                            '_PC', pc, '.rData'))
plot(envMap, main = substitute(paste('Maxent output, ', 
                                                        italic('Fraxinus quadrangulata '),
                                                        'occurrences')),
     sub = paste0('GCM = ', gcm, ', Buffer = ', buffer_distance))
plot(rangeMap, border = 'blue', add = TRUE)